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Abstract. 

The non-Gaussian distribution of primordial perturbations has the potential to 
reveal the physical processes at work in the very early Universe. Local models provide 

a well-defined class of non-Gaussian distributions that arise naturally from the non- 
linear evolution of density perturbations on super-Hubble scales starting from Gaussian 
field fluctuations during inflation. I describe the SN formalism used to calculate the 
primordial density perturbation on large scales and then review several models for the 
origin of local primordial non-Gaussianity, including the cuvaton, modulated reheating 
and ekpyrotic scenarios. I include an appendix with a table of sign conventions used 
in specific papers. 



1. Introduction 

The common presumption that primordial density perturbations have a Gaussian 
distribution is a powerful simplifying assumption that allows one to specify all the 
properties of the distribution once the two-point correlation function is known in real 
space, or equivalently the power spectrum in Fourier space. In particular the three- 
point and connected higher moments of the distribution vanish. On the other hand the 
statement that a distribution is non-Gaussian opens up an infinite array of possibilities. 
This has led to an assortment of empirical tests for non-Gaussianity of the primordial 
perturbations. By contrast there are relatively few non-Gaussian distributions that are 
motivated by theoretical models for the origin of structure in the early universe. While 
one can argue that a Gaussian distribution could describe density perturbations arising 
from a wide range of possible sources, any detection of deviations from a Gaussian 
distribution predicted by a specific theoretical model would be strong evidence in 
support of that model. 

Vacuum fluctuations in light, weakly-coupled scalar fields during a period of 
inflation (defined here as accelerated expansion, a > 0) in the very early universe provide 
a natural origin for an almost Gaussian distribution of field perturbations on large scales. 
Fluctuations in a free quantum field on small scales with comoving wavenumber, k, 
are swept up to scales much larger than the comoving Hubble scale, H~^/a = 1/a, 
which shrinks during infiation. On super-Hubble scales {k < aH) damping drives the 
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fluctuations into a squeezed state and they can effectively be treated as a classical 
random distribution Models of slow-roll inflation in the early universe are driven 
by canonical fields whose self-interactions are small relative to the Hubble scale and 
hence their distribution on scales close to the Hubble scale remains well described by a 
Gaussian distribution with vanishing third- and higher-order connected functions |21 13] . 

As inflation continues to stretch the perturbations up to scales far larger than the 
Hubble scale, spatial gradients are expected to become negligible and the we may treat 
the evolution locally as effectively "separate universes" [H E] whose initial conditions 
are set by the Gaussian distribution of the scalar field values during inflation. This 
motivates the study of local models of primordial non- Gaussianity where the distribution 
of primordial density perturbations, (, can be described by a local function of one or 
more Gaussian random fields, ({6(p'). A linear function of Gaussian fields is itself 
Gaussian, so a non-Gaussian distribution implies non-linearity and the study of non- 
linear cosmological perturbations. This local model for non-Gaussianity turns out to 
be a very good description of non-Gaussianity in some simple physical models for the 
origin of structure in the very early universe. 

Conventionally the primordial perturbations are characterised by the metric 
potential in the matter-dominated era, $ = (3/5)C- In Fourier space we define the 
power spectrum and bispectrum as 

($k,$k.) = {2nfP^{h)6''{k^ + ka) , (1) 
(^k.^k.^-ka) = {27rfB^{h, h)6%k^ + k2 + kg) , (2) 

and the amplitude of the bispectrum relative to the power spectrum is then given by 
the dimensionless parameter 

r (U u u \ = B<s,{ki, k2, h) 

^"""^^ ''~2 mh)p^{h) + p^{h)p^{h) + p^{h)p^{h)] • ^ ^ 

In a perturbative expansion, the second-order expression for the metric potential 
as a local function of a single Gaussian field, 0, is given by 

$(x) = 0(x) + /nl (0'(x) - (02)) + . . . . (4) 

In this special case /nl is, by construction, a constant parameter independent of spatial 
position or scale. However in more general local models the ratio of the bispectrum to 
the power spectrum given by Eq. ^ may be taken as a definition of /nl which is then 
scale and shape-dependent. Note that following Komatsu and Spergel [6] we adopt a 
sign convention for the metric potential $, and hence /nl, which is the opposite of that 
used by, for example, Mukhanov et al [7J and Maldacena j2j. Specific sign conventions 
are summarised in a table in an Appendix. 

In section 2 I review the 5A^-formalism that is commonly used to calculate 
the primordial density perturbation on large scales and its higher-order correlations, 
including the use of Feynman-type diagrams and how the 6N formalism extends to 
the description of primordial isocurvature density perturbations. In section 3 I review 
several examples of models for the origin of structure in the very early universe and 
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the local non-Gaussianity they give rise to. I conclude in section 4 with brief review of 
current observational constraints. 

2. The 5N expansion 

A powerful technique for calculating the non-linear primordial density perturbation 
in many cosmological models is the 6N formalism [HI El El |10], which identifies the 
primordial density perturbation with the perturbed logarthmic expansion, N = J Hdt. 

In Friedmann-Robertson- Walker (FRW) cosmology there is a preferred foliation of 
spatial hypersurfaces which are maximally symmetric, and on which the matter density 
and pressure are also homogeneous and isotropic. In an inhomogeneous universe we 
can use the difference between uniform-expansion hypersurfaces and uniform-matter 
hypersurfaces as a measure of inhomogeneity. At first order we define [HI [T2] 

Ci = Ci--V- (5) 
P 

where Ci and 5ip are the gauge-dependent spatial metric and density perturbations 
respectively [3 [12] (see Appendix A). This can be interpreted either as a curvature 
perturbation on uniform-density hypersurfaces (where 5p = 0) or a dimensionless density 
perturbation on uniform-curvature hypersurfaces (where C = 0). 

Local energy conservation equation [5] ensures that C^i remains constant for 
adiabatic density perturbations in the long-wavelength limit where (fc/aif)^$i — )■ 
and Ci approaches the comoving curvature perturbation, TZ. Thus during inflation it is 
convenient to evaluate the comoving curvature perturbation on large scales, k/aH — )■ 0, 
in order to determine the primordial density perturbation from single-field inflation 

Ci ^ 7^l = Ci - -6iif . (6) 



Because the non-adiabatic decaying mode rapidly decays outside the Hubble scale [I3] 
it is sufficient to evaluate the primordial perturbation due to scalar field perturbations 
shortly after Hubble exit, k = aH. However in multi-field inflation, additional light 
fields can lead to non-adiabatic pressure perturbations on super-Hubble scales and Ci 
becomes time-dependent in general [TH [T5] . 

Beyond linear order, we can identify the non-linear metric perturbation which 
coincides with (i at first order and on large scales [16] (see also [171 [18] and the article 
by Langlois and Vernizzi elsewhere in this issue [19] ) 

1 /-pC*.") dp 

hit) p + P 

where p{t) is the homogenous background density, p{t, x) denotes the local 
inhomogeneous density and the local expansion is given by iV = / Hdt, where H{t, x) = 
is the local Hubble expansion rate along comoving worldlines, = dx^/dr. 
A local form for the distribution of C, naturally arises from the evolution of matter 
fields on scales much larger than the Hubble scale in the very early universe. After 



1 rpit,^) rj n 

C(t,x) = 5Ar(t,x) + -/ (7) 
o Jm + P 
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averaging on some scale, L ^ H~^, such that spatial gradients and anisotropy can be 
neglected, then the local expansion is well-described by the Friedmann equation for a 
homogeneous universe with the corresponding local density and pressure. This is known 
as the "separate universe" approach [5] , and can be derived from the full inhomogeneous 
equations of motion as a long- wavelength limit in a gradient expansion [H [TTl dB] . In 
this approach the local evolution, along a given worldline, is determined by the initial 
values of the matter fields on an initial spatial hypersurface. 

Therefore one can evaluate ({tf,x), defined by Eq. ([T]), as 5N, the perturbed 
expansion on a uniform-density hypersurface (p(tj,x) = p(t/) for all x) at some final 
time tf after inflation has ended by evaluating the integrated expansion from some 
initial spatially flat hypersurface, 6N(ti,x.) = for all x [SI El [10]. In the separate 
universe framework [3] the inhomogeneous universe is modelled as a patchwork of locally 
homogeneous regions which, due to causality, evolve independently. In particular one 
can determine the non-linear local expansion N using the Friedmann equation for the 
local expansion as a function of the initial local flelds, ip\ti,x.) = (p{ti) + 6ip{ti,x.). 
Assuming N{(p^) is an analytic function of the initial field values, we have the Taylor 
series expansion [TU] 



where Nj = dN/d(p', etc, and summation is implied over repeated indices. If the 
field perturbations themselves are constructed as an expansion in a small perturbation 
parameter 



then we can write, order by order in a perturbative expansion = (^i + C2/2 + Cs/S . . ., 
where 



To calculate the statistics of the primordial perturbation, we therefore need to 
know the dependence of the large scale expansion upon the initial field values, N{ip^), 
and the statistical distribution of the initial field values. The distribution of the fields 
is determined by the quantum fiuctuations of the vacuum in any given infiationary 
model. We take the first-order perturbations to describe the vacuum fiuctuations of 
the non-interacting free fields, and build up the higher-order terms from the interaction 
Hamiltonian, following Maldacena |2] (see also the article by Koyama elsewhere in this 
issue [20] )• 

In canonical slow-roll infiation the linear field perturbations defined on unperturbed 
spatially fiat hypersurfaces [2Tl|22] are given by the fiat spacetime vacuum on sub-Hubble 
scales and are then approximately constant on super-Hubble scales with power spectra 



C = iV(/) -N = Nj5ip' + -Nij5ip'5ip' + -NuKSif'dip'dip'' + . . . , (8) 

2 D 




(9) 



Ci 

C2 



(10) 

(11) 
(12) 



(13) 
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where P{ki) ^ Hl/2k\ and if* is the Hubble scale at Hubble-exit, k = a^H^,. Note that 
the dimensionless power spectrum, related to the variance in real space, is multiplied 
by a volume factor in fc-space to give 

nk) = i^Mk) ^ (T] . (14) 



<1. (15) 



(27r)3 ' ' V27r, 
and this is approximately scale invariant, with spectral tilt 

dlnV 
dink 

Scalar field interactions are suppressed by slow-roll parameters and so field 
perturbations on scales close to the Hubble scale are approximately Gaussian. In the 
following we will assume that all connected higher-order moments of the scalar field 
perturbations vanish at some fixed time during inflation, shortly after all relevant scales 
have crossed outside the Hubble scale, k < aH. Hence N is well-described by a local 
function of Gaussian random fields soon after the smallest relevant modes have left the 
Hubble scale [fj 

We will assume that S2(p^ , S^ip^ , etc, vanish sufficiently close to Hubble exit and 
thus the higher-order moments are given by the first terms on the right-hand-side of 
Eqs. (11) and (12). In particular the primordial power spectrum is given, at leading 
order, by 

(Ck.Ck.) = NjNj{6^i^6^Q = {2n)-'P^{k,)6'{k^ + k^) , (16) 

where 

P^ik,) = NjN'5''P{k,) , (17) 

In the next subsection we will discuss higher-order contributions to the primordial power 
spectrum coming from higher-order terms in the 6N expansion. These can be identified 
as "loop" corrections in a diagrammatic approach where the leading order result 
(16) is the "tree- level" term. 

Similarly we can construct the leading order terms in the bispectrum and 
trispectrum 

(Ck.Ck.Cks) = {2nfB^{k,, k2, ks)6^{k, + k2 + ks) , (18) 
(Ck,Ck.Ck3Ck4) = (27r)3Tc(fci, k2, ks)5'ik^ + ka + kg + k4) , (19) 

where 

Bdku k2, ks) = NiNjN'' [P{ki)P{k2) + 2 perms] , (20) 
Tdki, k2, k,) = NijkN'N'N'' [P{k2)P{k,)P{k^) + 3 perms] 

+ NijN'^N'N'' [P{ks)P{k^)P{k,s) + 11 perms] , (21) 

and ki2 = |ki — k2|, etc. 

I Alternatives to slow-roll models of inflation based on non-canonical scalar field Lagrangians with non- 
linear kinetic terms, such as k- inflation [23] or DBI inflation 24J, can lead to significant non-Gaussianity 
of the field perturbations at Hubble exit and hence produce different types of non-Gaussianity, different 
from the local form. 
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(a) (b) 



Figure 1. Power spectrum diagrams for SN expansion witli Gaussian fields (a) at tree 



level, and (b) one loop, corresponding to each term in Eq. (22 1 



2.1. Diagrammatic approach 

The perturbative expression of the non-hnear expansion as a function of the field values 
during inflation ([s]) is conveniently expressed in terms of Feynman-type diagrams |2S] 
(see also [25]). Similar diagrams have been used in the analysis of large-scale structure in 
cosmology [23 12H] • Every term in the n-point function can be identified with a diagram 
with n external lines. To perturbative order r one should draw all diagrams with r 
internal propagators. Tree-level diagrams correspond to r = — 1 while diagrams with 
r > n are loop corrections. If the expansion is a local function of Gaussian fields then 
the internal propagators have no self-interactions. Examples are shown in Figure 1. A 
comprehensive prescription for drawing diagrams and constructing the corresponding 
terms to any given order are given by Byrnes et al [26] . 

Here we will consider only diagrams for Gaussian fields, ip^. The power spectrum 
up to one-loop is then shown in figure 1, corresponding to 

^ J d'qd'q' (^^NABNcDS^q! - q + k) + NANBcnS'iq! - k)) 
5^^(5^^P(g)P(g') . (22) 



+ 



Similarly the bispectrum and higher-order correlations can be evaluated including an 
arbitrary number of loop corrections. 

Because the loop corrections include integrations over all internal momenta, and 
the scalar field perturbations are expected to have almost scale invariant power spectra 
during inflation, the loop corrections can give corrections that, although only slowly 
growing on large scales, formally diverge in the infra-red limit. This issue is discussed 
by Seery elsewhere in this issue [29]. Typically the corrections remain small in a region 
which is not much larger than the observable universe, but in some cases, for instance if 
the tree-level terms are absent [3U], then loop corrections may play an important role. 



Some loop corrections, such as the second one-loop term in Eq. (22), can be 



identified as dressing the vertices. Byrnes et al [22] noted that these can be interpreted 
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as renormalising the derivatives of N in terms of the average local derivatives: 

(Na) = J d\NA + ^ j d'qNABcS''''P{q) + • ■ ■ (23) 

Renormalising the vertices in this way simplifies the expressions for the loop corrections, 
to yield only one new term at each new order. In real space this is sufficient to resum all 
divergences, but in Fourier space the dependence on IR wavenumbers on the convolution 
remains. For instance the power spectrum to one loop becomes 

+ (2^ / d'cil{NAB){Ncn)6^^S^''P{k)P{k - q) . (24) 
2.2. Single-field vs multi-variate local models 

If the primordial density perturbation ([8| is a local function of a single scalar field 



C = + \n"5^^ + In"'5^'' + . . . , (25) 
2 o 

then the expressions for the higher order moments of the primordial distribution, e.g.. 



(20) and (21), simplify considerably and we can write 



Bciki, k2, ks) = ^/nl [Pdki)Pdk2) + 2 perms] , (26) 
54 

(fci, k2, ks) = ^Qnl [Pc(/c2)^c(^3)^c(^4) + 3 perms] 

+ Tnl [Pdk3)Pdk4)Pdki3) + 11 perms] , (27) 
where at leading order (tree-level) the non-linearity parameters are given by [TOl l3T] 
5 A^" 36 25 A^'" 



The 6N expression for the primordial perturbation ([8j) goes beyond the simplest 
local model ^ by considering the case where the primordial perturbation is a function 
of more than one Gaussian random field. Including terms up to second order we have 

$(x) = ^ <l>\x) + Y: fiJ ($'(x)$^(x) - ($^<l>^)) + . . . , (29) 

where we have defined $^ = {3/5)Nj6ip' and fu = {6/5)Njj/NjNj (no sums). The 
non-linearity parameter /nl defined by Eq. ^ is then fc-dependent and given by 



r J:I,J Njj [Pi{kl)Pj{k2) + Pl{k2)Pj{k,) + Pi{k3)Pj{k,)] 

N,Nj [Pj{k,)Pj{k2) + Pi{k2)PAk,) + Pj{k,)Pj{k,)] ■ ^ > 

This is a constant if the scale-dependence of all the fields which contribute to the 
primordial perturbation is the same, but in general it is scale and shape dependent, 
/c-dependence of /nl defined by Eq. ^ also arises due to interactions and hence 
non- Gaussianity of fields on super Hubble scales [3^, but this goes beyond the strict 
definition of local models considered here. If the fields are all approximately massless 
during infiation then their scale-dependence and interactions are small and the k- 
dependence of /nl is expected to be weak (but see 
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2.3. Isocurvature non- Gaussianity 

Having identified the non-linear primordial density perturbation, (, with the perturbed 
expansion, 6N, up to a uniform density hypersurface, it is straightforward to extend 
this to n primordial density perturbations in an n-component system, where we define 

O J Pa Pa ' J^a 

Relative density perturbations between the different components which leave the 
total energy density unperturbed are know as isocurvature perturbations. In particular 
a relative perturbation between photons and baryons around the time of last scattering 
of the CMB photons corresponds to a perturbation of the local entropy (determined by 
the photon number density) per baryon and hence such perturbations are also known 
as entropy perturbations. At linear order we have 

^B = ^-^. (32) 

We can identify this with the difference between the two primordial density 
perturbations (b and C7, up to a factor of 3. Hence beyond linear order we define 



Sb = 3 (C^ - Cb) . (33) 

A similar definition is used for the cold dark matter isocurvature perturbation, Scdm, 
and the neutrino isocurvature perturbation, 5*,^, relative to the photon number density. 

More generally for an n component system we can decompose an arbitrary density 
perturbation into one adiabatic density perturbation, C = Ca for all components a, and 
n — 1 independent isocurvature modes 

Sap = 3(C - C/3) • (34) 

One can define non-linearity parameters for the isocurvature bispectrum at leading 
order analogous to the adiabatic bispectrum ^ 

r , , N Bs{ki,k2,h) 

fs[l^„ K,) 2 [p^^k,)Ps{h) + Ps{h)Ps{h) + Ps{ks)Psih)] ■ ^^'^ 

If the isocurvature modes are uncorrelated with the adiabatic density perturbation, and 
we use the linear relation between primordial matter isocurvature perturbations and the 
Newtonian potential, $ = Scdm/5, we obtain 



r ^ 5/5 Psik^)Psik2) + Psik2)Psik,) + Psik,)Psik,) 
^""^ 162 P^ik^)P^ik2) + P^{k2)P^{k^) + P^{k^)P^{k^) ■ ^ ' 

However one should be wary of interpreting bounds on the adiabatic /nl as any 
constraint on isocurvature non- Gaussianity which leads to a distinctive bispectrum on 
the CMB sky [3H |36] . Moreover the isocurvature modes are in general correlated with 
adiabatic modes [37] and hence there may exist non-vanishing cross-correlations, {S(() 
and (^^C) [35]. 
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3. Models of non-Gaussianity 

3.1. Single inflaton 

A homogeneous scalar field cosmology evolves like that of a fluid with time dependent 
density and pressure, but an inhomogeneous scalar field is not in general a barotropic 
fluid [38] • However the non-adiabatic part of the field's pressure perturbation is 
proportional to the field's comoving density perturbation [13]. The Einstein equations 
energy and momentum constraints require that the comoving density perturbation 
vanishes in the large scale limit where the Newtonian potential stays finite [12] and 
hence if the local energy density and pressure during inflation is dominated by a single 
scalar field then perturbations must be adiabatic in this large scale limit. The field 
perturbations enter a squeezed state in phase space and the local time derivative of the 
field is no longer independent of the local field value. 

This makes it straightforward to predict the primordial density perturbation and 
hence the local non-Gaussianity for inflation models with a single inflaton field. As 
all perturbations are adiabatic then the local density and pressure follow the same 
phase-space trajectory in phase-space as the background cosmology with 5p/5p = p/p. 
This ensures that the perturbation ( is non-linearly conserved in this large scale limit 
[SI EHl 113 [ISl [IE] and we can calculate the primordial density perturbation, long after 
inflation, in terms of quantities during inflation. Using the 5A^-formalism we just need 
to calculate the effect of field perturbations on the local expansion, N, during inflation. 

To first order we have 

TT 

Ci = N'6<^ = 6<^ , (37) 

where 6ip denotes field perturbations in the spatially flat gauge. This is time-independent 
on large scales, but in practice it is convenient to evaluate it soon after Hubble exit 
[k = aH) where the amplitude of the field perturbations are given by P{k) ~ H^/2k^. 
We then obtain the simple result for the primordial power spectrum at leading order 

Pdk) = (I) Pik) . (38) 
At second order, assuming a Gaussian field 6(f, we have 

G = A"V=(-|^ + ff)V. (39) 

In terms of the usual slow-roll parameters e = —H/H"^ and rj = —{(p/H(p) — [H / H"^) 
we have 

G = (2e - r/) Cl ■ (40) 

Thus we see that the second order primordial density perturbation is suppressed with 
respect to the square of the first-order perturbations and hence |/nl| ^ 1 [lOlIll]. In fact 
at this order in the slow-roll approximaton we can no longer neglect the intrinsic non- 
Gaussianity of the fields at Hubble exit [2], and a simple argument where we consider 
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long-wavelength perturbations to set the local Hubble expansion as shorter wavelengths 
leave the Hubble scale ^] shows that in the squeezed limit, where kf <^ A;| + we 
find 

hLiki,h,h) = l{3e-r]) . (41) 
o 

More generally, for any single field model producing adiabatic perturbations during 
infiation the bispectrum can be related to the scale-dependence of the power spectrum 
in the squeezed limit and we have [13] 



hL{kuk2,h) = ^{l-n^). (42) 

where — 1 = cilnP^/ciln /c. Observational bounds on the spectral tilt, \n^ — 1| <^ 1, 
thus ensure that local-type non-Gaussianity is very small in single-infiaton field models. 

3. 2. Curvaton 

The curvaton scenario [HI HHl HHl El SH] provides a simple model for the origin of the 
primordial perturbations which could exhibit significant non-Gaussianity of the local 
type m- 

The curvaton is a weakly-coupled scalar field which is light during inflation, m H, 
but it has a negligible energy density at that time. Quantum fluctuations of the fleld 
during inflation generate an almost scale invariant spectrum of fleld perturbations on 
super-Hubble scales. We assume that the weakly-coupled fleld remains decoupled from 
the inflaton and its decay products at the end of inflation. The fleld begins to oscillate 
about the minimum of its potential when the Hubble rate drops below the mass of the 
fleld some time after inflation. An oscillating massive fleld, m ^ H, has a pressureless 
equation of state (averaged over several oscillation times) and hence the curvaton energy 
density decays as a~^, but grows relative to the energy density of radiation, oc a~^. This 
is the Polonyi or moduli problem of weakly coupled scalar flelds which can come to 
dominate the energy density of the early universe, disrupting the conventional hot big 
bang. This is not a problem so long as the moduli decay and their decay products 
thermalise before the epoch of primordial nucleosynthesis. The inhomogeneous density 
of the curvaton is transfered to the radiation when the curvaton decays leading to a 
primordial density perturbation on super-Hubble scales. 

Note that, in contrast to inflaton perturbations, the curvaton fleld fluctuations are 
isocurvature fleld perturbations during inflation and thus give rise to non-adiabatic 
pressure perturbations after inflation which leads to a change in the value of the 
perturbation, (, on super-Hubble scales. Thus the bispectrum in the squeezed limit 



is not constrained by the relation (42) for perturbations generated by the adiabatic 
mode during inflation. 

In the simplest example of a curvaton fleld, with potential V{x) = 
the curvaton fleld is a free fleld during inflation with no self-interactions and hence is 
described by a Gaussian fleld. The quadratic potential naturally leads to a primordial 
density fleld that is a quadratic local function of the Gaussian curvaton fleld and 
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hence the simplest curvaton model is well described by the simplest local model of 
non- Gaussianity given in Eq. Q. 

The density perturbation for the oscillating curvaton field on spatially flat 



hypersurfaces is given by Eq. (31) 



3 Px 

In the case where the curvaton density remains sub-dominant throughout we can 
assume a linear transfer to the primordial density perturbation, with transfer efficiency 
r ~ ^^x.decay = Px /Ptot | decay , and we can write 

C = rCx = Ci + ^/nl (Ci' - (Ci')) . (44) 
where we identify 

= (45) 

3 X 



and the non-linearity parameter 

h. = l. (46) 

For r <^ 1 we can have large, positive /nl|§} In this case current observational bounds 
on /nl 163] place a lower bound on r > 0.01. 

For r ~ 1 we need to include the gravitational effect of the curvaton density when 
it decays. If we assume an instantaneous decay hypersurface given by the local Hubble 
rate = F, we obtain the non-linear relation between the curvaton perturbation 
and the primordial perturbation ( 

(1 ^Xidccay)^^ ^Xi<iscayC ^'^ 1 , (47) 



from which we recover Eq. (45) at hnear order, where r = 3f2^^decay/(4 — ^^x.decay), and 
at second and third order we obtain the non-linearity parameters as a function of r 

/N. = ^(l-y-^j. (48) 
*c^-?fl-i-^-^V (49) 



Qr \ 18 9 3 ^ 

A measurement of both /nl and (7nl would be an important consistency test of the 
simplest curvaton model as both are a function of a single parameter, r. For a sub- 



dominant curvaton, r ^ 1, we recover Eq. (46) for /nl and (7nl — — 10/nl/3 [51j. In 



§ There has been some confusion in the hterature over the "correct" sign for /nl in the curvaton 
scenario. See Appendix A for a table of sign conventions used in specific papers. In my own papers, 
the original result for /nl in Ref. [49] had the "correct" sign (that used here, which coincides with that 
used for observational constraints ^) due to a combination of the sign convention used and a sign error 
in the arithmetic. The mistake in the arithmetic was corrected, for instance in Ref. [50j . but then the 
sign was "wrong" because of the sign convention used. 
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the opposite limit, as r — )■ 1, where the curvaton dominates the energy density of the 
universe before it decays, we have /nl — )■ — 5/4 and (7nl 25/12. 

Other possibihties that have been studied include a kinetic-energy dominated era 
(rather than radiation dominated) before curvaton decay [53] and multiple curvaton-type 
fields Even in the case of multiple curvaton fields contributing to the primordial 
density perturbation we recover the same lower bound /nl > —5/4 as in the single 
field case [SHI ES]- A much wider range of local non- Gaussianity becomes possible if 
we consider the effect of self-interaction terms in the curvaton potential, such that 
V = m?x^ /2 + Ax" + • • • with n > 3 which can lead to non- linear evolution of the field 
before it decays [STl ESI EH ED]- In this case one may find large negative values of /nl 
for specific model parameters and initial conditions. 

Another distinctive aspect of the curvaton model is the possibility to leave residual 
isocurvature perturbations after the curvaton has decayed |1H HSl HT] . Non-adiabatic 
perturbations on super-Hubble scales during inflation are necessary, but not sufficient, 
condition for the existence of primordial isocurvature perturbations after inflation. 
Whether or not residual isocurvature modes survive after the curvaton decays depends 
upon the reheating history [HI]. If all the curvaton decay products thermalise with 
vanishing chemical potential then the primordial perturbation must be adiabatic [62] . 
But if a baryon or lepton asymmetry has already been created, or the curvaton decay 
itself produces the asymmetry, or the dark matter has already decoupled, then a residual 
isocurvature perturbation, Sx-, may be left, and it will be completely correlated (or anti- 
correlated) with the total density perturbation, C,. 

For example, if the curvaton decay is the out-of-equilibrium process which breaks 
time-reversal invariance and violates the baryon number then the baryon number inherits 
the same density perturbation as the curvaton, = Cx? orders [35]. In this case 

bounds on the amplitude of the linear primordial isocurvature perturbation [03] requires 
that the curvaton dominate the energy density when it decays, r ~ 1, so that ~ 
and, from Eq. (33), 5^ ^ C-B- In this case the intrinsic non-linearity parameter of the 
isocurvature perturbation is of order unity, and the contribution to the CMB bispectrum 
is suppressed relative to that from the total density perturbation. 

A mixed curvaton-infiaton scenario yields a much richer phenomenology of 
adiabatic, non-adiabatic and non-Gaussianity of primordial perturbations [Ml EH] EEJ 
[6711351 EH. 



3.3. Multiple-field inflation 

The curvaton is just one example of how fields other than a single infiaton driving 
inflation, could play a significant role in determining the primordial density perturbation 
after the end of inflation. The presence of non-adiabatic perturbations on super-Hubble 
scales allows in principle for the evolution of the large-scale density perturbation after 
Hubble-exit [12] and hence local non-Gaussianity in the primordial density perturbation 
some time after inflation. 
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One of the best studied models is double inflation with two massive but non- 
interacting fields [8[ |37] where the non-linearity parameter though not necessarily slow- 
roll suppressed is not expected to be much larger than unity [691 U0\ 171] . Primordial 
isocurvature perturbations have also been studied in this model where one of the massive 
fields is identified as a dark matter candidate [721 EZl ES] • 

In general it seems surprisingly difficult to produce large non-linearity parameters 
during slow-roll inflation, but examples can be constructed ^J. For more discussion 
see the article by Tanaka et al in this volume [71] , and the recent review by Byrnes and 
Choi [75]. Non-Gaussianities could be produced by the breakdown of slow-roll due to 
features in the potential [75] or particle production [77]. 

In inflationary models where more than one field remains light until the end of 
inflation it is necessary to consider possible evolution of the density perturbation on 
large scales at or after the end of inflation. The end of slow-roll inflation is an epoch 
at which the large-scale density perturbation may have a non-linear dependence upon 
non-adiabatic modes, especially it there is an abrupt change in the equation of state 
[781 [7911801 EI]. 

It is very natural to realise this possibility in hybrid inflation models [H2] where 
inflation is ended by a tachyonic instability triggered in a waterfall field, leading to 
a rapid phase transition and decay of the false vacuum. In string theory models of 
inflation this may describe the collision of branes in a higher-dimensional space [83] . 
The instability corresponds to a surface in field space, and if the inflaton trajectory is 
not orthogonal to this surface then non-adiabatic field perturbations, orthogonal to the 
trajectory [13] , are converted into density perturbations at the phase transition. In the 
simplest case, taking the extreme slow-roll limit and assuming instantaneous reheating, 
the energy density will be almost uniform on the spatial hypersurface corresponding 
to the phase transition, and thus the perturbed expansion to this surface, 6N, can be 
identified with the primordial density perturbation, ( in Eq.([7]). 

3.4- Modulated decay 

A specific example where the primordial density perturbations may be produced just 
after the end of inflation is when the decay rate of the inflaton is a function of one or 
more moduli fields [SH |85]. Because this occurs some time after Hubble exit during 
inflation, this produces local type non-Gaussianity [86l [871 EH] . 

If we approximate the inflaton reheating by a sudden decay, as we did earlier for 
the curvaton decay, we find an analytic estimate of the non-linear density perturbation. 
In the case of modulated reheating, the decay occurs on a spatial hypersurface with 
variable local decay rate and hence local Hubble rate H^i = r(x). Before inflaton 
decay, the oscillating field has a pressureless equation of state and there is no density 
perturbation perturbation. Setting = in Eq. ([T]) thus gives the perturbed expansion 
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on the decay hypersurface 

6N, = - 




(50) 



After the decay the we have radiation with equation of state p = p/3 and hence from 
Eq. ([7]) a density perturbation 



(51) 



Ehminating SN^ and using the local Friedmann equation, p oc H^, to determine the 
local density in terms of the local decay rate, Hd = r(x), we have 




c 



At first order we recover the linear relation 



At higher order we obtain 



C2 




1^ 



and hence 



/i 



NL 



5 1- 




r'T\ 



(52) 

(53) 

(54) 
(55) 

(56) 
(57) 



and so on to higher order. Thus for an approximately linear modulation function T{x) 
we have /nl — 5 and (/nl — 100/3 while for F oc we obtain /nl — 5/2 and (^nl — 25/3. 
Similar results are found if the mass of the infiaton decay products is modulated, so- 
called inhomogeneous mass domination [891 EZ] • 

Resonant decay or preheating may be sensitive to non-adiabatic modes [901 El 1221 
[93I IHH IHH]- In these models observational constraints on primordial non-Gaussianity 
impose significant constraints on the allowed parameter values and initial conditions 
For a review of infiaton dynamics and reheating see Bassett et al 



3.5. Ekpyrotic model 

There have been numerous attempts to construct alternatives to inflation as a model for 
the origin of primordial perturbations. Many of these exploit the similarities between an 
inflationary expansion and a collapse phase where the comoving Hubble scale decreases 
and thus quantum vacuum fluctuations evolve into the super-Hubble regime [^Sl EZ]- 
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As the Hubble rate grows during collapse, these models generally require some form of 
instability for the perturbations to grow on super-Hubble scales and maintain a scale- 
invariant spectrum. For example, a simple collapse model driven by a scalar field with 
pressureless equation of state generates a scale-invariant spectrum for ( on super-Hubble 
scales [23 EB]- The non-Gaussianity in this has recently been calculated [HS] and shown 
to yield /nl ~ 1, though not of the local form in Eq. Q. Ultimately the four-dimensional 
low-energy effective theory must break down as the collapse rate approaches the Planck 
scale and this leaves some uncertainty about how the perturbations are transfered to 
the expanding hot big bang, but at least within the collapse phase it is possible to study 
the growth of perturbation ( and its non-Gaussianity. 

Recent attention has focussed on the ekpyrotic model |100^ 1101] . In the 4D effective 
theory this corresponds to a cosmology driven by scalar fields with steep exponential 
potentials, Vi{(p^) cx exp(— c/y?'^) (no sum) leading to an ultra-stiff equation of state 
p/p ^ 1. The separate universe picture is an excellent approximation on super- 
Hubble scales during an ekpyrotic phase [391 HHj - Thus, for fields with canonical 
kinetic Lagrangians, the local values of the fields perturbations at Hubble-exit have 
an approximately Gaussian distribution and set the initial conditions the subsequent 
local expansion history. As a result we expect the primordial non-Gaussianity to have 
a local form, and due to the strong self-interaction terms in steep potentials, we expect 
the non-linearities to become large |103] . 

The adiabatic mode during an ekpyrotic collapse leads to a steep blue spectrum 
for ( |104] but in the presence of two or more fields with steep exponential 
potentials it is possible to produce an almost scale-invariant spectrum of isocurvature 
perturbations |1U5[ I1U6[ I1U2[ I1U7[ I1U3] . Consider the simple case of two fields with 
potential 

V = -Vi exp (-ciy^^i)) - V2 exp (-C2</^('^) , (58) 
Performing a rotation in field-space this can be re-written as |lU8j 



V = Vq exp(— CO") 



2 -^sc^ 3 K^c^ 4 

-1 s i=S s 

2 12v^ 96 



(59) 



where the fast-roll parameter = c^^ + C2 ^ and we identify s as the tachyonic field 
direction which acquires an almost scale-invariant spectrum of perturbations about s = 
on super-Hubble scales, with spectral tilt 

d\nV 4 , , 

For an almost scale- invariant spectrum we require 3> 1. 

Note that in the 6N formula the scale dependence of the primordial power 
spectrum (at leading order) follows directly from the scale-dependence of the initial 
field perturbations, — 1 = 4/c^, but the amplitude depends on the subsequent 
expansion history. In the ekpyrotic scenario different models have been proposed for 
the conversion of the isocurvature field perturbations to density perturbation (. The 
tachyonic instability itself will lead to a phase transition to an ekpyrotic phase dominated 
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by just one of the fields, |109j . In this case N ((p^^\ Lp^'^^] 
for the simple model (58) and we find |110t llllj . 



/i 



NL 



5 2 25 4 



can be calculated analytically 



(61) 



12 ' ' " 108 

Note that c| > 3> 1 and therefore the non-linearity parameters are expected to be 
large. Tight observational bounds on a negative value for /nl > —10 [63] require cj < 24 
in this simple model and hence from Eq. (60) an unacceptably blue primordial spectral 
tilt n(^ — l> 0.16. An acceptable tilt and non-Gaussianity would require modifications to 
the simple potential Eq. (60) which itself could trigger a transition and a isocurvature- 
curvature conversion. The amplitude and sign of the non-linearity parameters then 



becomes model-dependent but one generally finds 



NL 



IS'nlI 



Alternative conversion mechanisms would produce a different amplitude of 
primordial perturbations and hence non-Gaussianity. In models where the conversion 
happens at a reflection in field space during a kinetic dominated phase give estimates 
for the non-linearity parameters |111] . 



/i 



NL 



2^2 



K3C + 5, 



5'NL 



-20c^ 



32 



24 



(62) 



K3 and are expected to be of order unity. Again the non-linearity parameters are 



expected to be large with |(7nl 
such that I/nl 



/nl ~ c^. In this case we see that even if is small, 
1 [corresponding to — ^ 1/c in Eq. (58)], then we would expect 
Ig^hl to be large, unless is also unexpectedly small. We note the further possibility 
that isocurvature field fluctuations are converted to density perturbations only after the 
bounce |113j . as in curvaton or modulated decay models. 



4. Conclusions 



The distribution of primordial density perturbations contains more information than just 
the power-spectrum. Among a plethora of different possible forms of non-Gaussianity 
the local models provide a clearly defined class which arise naturally from the evolution 
of the density perturbation on super-Hubble scales from initially Gaussian vacuum 
fluctuations during inflation. In particular the (2 + n)-point function of a non-Gaussian 
distribution generated by a local function of a single Gaussian field can be described at 
leading order by the two-point function and n non-linearity parameters, /nl, fi'NL, etc. 

Non-Gaussianities have the potential to reveal the physical interactions at work in 
the very early Universe. For adiabatic perturbations during inflation the non-linearity 
parameters are related to the scale-dependence of the power spectrum [2l |13] and hence 
must be small. Thus any detection of local non-Gaussianity of the primordial density 
perturbation would be evidence of non-adiabatic perturbations on super-Hubble scales, 
e.g., the presence of multiple light fields during slow- roll inflation. 

Current observational limits on /nl in the single-field local model come from 
combining WMAP cosmic microwave background bounds (—10 < /nl < 74 [63]) with 
large-scale structure (—29 < /nl < 70 |114j ). yielding —5 < /nl < 70 at 95% confidence 
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limit [63]. Recent work also places a bound —3.80 < S'nl/IO^ < 3.88 from the WMAP 
5-year data |115] . Care should be taken when applying these bounds to non-local or 
even multi-variate local models. CMB bounds are obtained using matched filtering 
techniques to construct an optimal estimator for /nl defined by Eq. Q. It estimates 
the amplitude of a second-order primordial perturbation which has a specified correlation 
with the first-order perturbation. On the other hand bounds from large scale structure 
arise from the form of the (general relativistic) Poisson equation |116] which implies 
that a local form for the primordial Newtonian potential requires a non-local form for 
the comoving density contrast, leading to a distinctive scale-dependent bias |117j . 

There is great potential for future discovery coming from future CMB missions 
such as ESA's Planck sateUite which has the potential to bound |/nl| < 5 and large- 
scale structure surveys, including 21cm line radio surveys |118] . There is also plenty 
of scope for further theoretical developments given the data that already exists to test 
non-Gaussian models. New theoretical templates will need to be derived for optimised 
constraints to be placed on alternative theoretical models of non-Gaussianity. At some 
point primordial non-Gaussianity will be discovered as it inevitably arises from the non- 
linear evolution of density perturbations and the last-scattering of the CMB photons 
[119[ 1120] , and it will then be a question of using all the available data to disentangle 
the different contributions. 
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Appendix A. Sign conventions 

We collect together here definitions (at linear order) and sign conventions for the gauge- 
invariant curvature perturbations used in the literature. 

Given the FRW line element with scalar metric perturbations 

ds'^ = a'^iri) [-(1 + 2A)dt^ + 2diBdx'dri+ 

{{l + 2C)6ij + 2didjS}dx'dx^] , (A.l) 

the intrinsic spatial curvature is given by the Ricci scalar on constant-?] hypersurfaces 

i^)R=-±d^C. (A.2) 
We can define three gauge invariant curvature perturbations at first-order: 
(i) uniform-density curvature perturbation 



Local non- Gaussianity from inflation 



18 







7^ 




/ 1\TT 
J INlj 


this article 


+ 


+ 


+ 


+ 


Malik & Wands ^ 








+ 


Bassett, Tsujikawa & Wands |50] 










Lyth & Rodriguez |10j 


+ 








Lytli, Malik & Sasaki (sections 3 & 4 only) |16j 










Maldacena [2] 










Komatsu & Spergel [B] 






+ 


+ 


Liddle & Lyth [1:^1] 




+ 






Mukhanov, Feldman & Brandenberger [7j 




-c 






Kodama & Sasaki |122j 




+ 


+ 




Bardeen |123 




+0m 







Table Al. Summary of sign conventions and alternative notations used in specific 
papers. 



(ii) comoving-orthogonal curvature perturbation 

Tl^C-^, (A.4) 

(iii) longitudinal-gauge curvature perturbation 

<l> = C + 'H{B-£'), (A.5) 

where primes denote derivatives with respect to the conformal time t] and "H = a' /a. 
Note that in the longitudinal gauge the curvature perturbation can be simply related 
to the Newtonian potential. 

The sign of the metric perturbation at the gauge-invariant variables formed from 
it are arbitrary and different authors have adopted different signs and notations. A 
different choice of sign for $ in Eq.(|2| then leads to a different choice of sign for /nl in 
Eq.([3]). This too is purely conventional, but once this choice is specified observations 
are sensitive to the sign of the resulting $ and hence the sign of /nl- Using the sign 
conventions in this paper, a CMB sky with positive /nl has more cold spots, and one 
with negative /nl has more hot spots. 

Different notations and sign conventions used in specific papers are listed in 
table Al. 



References 

[1] D. Polarski and A. A. Starobinsky, Class. Quant. Grav. 13, 377 (f996) |arXiv:gr-qc /9504030| . 

[2] J. M. Maldacena, JHEP 0305 (2003) Of 3 ar Xiv:astro -ph/02 f 0603 ] . 

[3] D. Seery and J. E. Lidsey, JCAP 0509 (2005) Off farXiv:astro -ph/0506056| . 

[4] D. S. Salopek and J. R. Bond, Phys. Rev. D 42, 3936 (1990). 

[5] D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle, Phys. Rev. D 62, 043527 (2000) |arXiv:astro-" 
'ph/0003278|. 

[6] E. Komatsu and D. N. Spergel, Phys. Rev. D 63, 063002 (200f) |a rXiv:astro-ph/0005036] . 



Local non- Gaussianity from inflation 



19 



[7] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215, 203 (1992). 

[8] A. A. Starobinsky, JETP Lett. 42, 152 (1985) [Pisma Zh. Eksp. Teor. Fiz. 42, 124 (1985)]. 

[9] M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95 (1996) 71 |arXiv:astro-ph/9507 00l . 

[10] D. H. Lyth and Y. Rodriguez, Phys. Rev. Lett. 95 (2005) 121302 |arXiv:astro-ph/0504045|. 

[11] J. M. Bardeen, DOE/ER/40423-01-C8 Lectures given at 2nd Guo Shou-jing Summer School on 

Particle Physics and Cosmology, Nanjing, China, July 1988 (1988). 

[12] K. A. MaHk and D. Wands, Phys. Rept. 475, 1 (2009) arXiv:0809.4944 [astro-ph]]. 



[13] C. Gordon, D. Wands, B. A. Bassett and R. Maartens, Phys. Rev. D 63, 023506 (2001) |arXiv:astro- 



ph/0009131 



[14] A. A. Starobinsky and J. Yokoyama, arXiv:gr-qc/9502002 



[15] J. Garcia-BeUido and D. Wands, Phys. Rev. D 53, 5437 (1996) [arXiv:as tro-ph/9511 029]. 

[16] D. H. Lyth, K. A. MaHk and M. Sasaki, JCAP 0505 (2005) 004 arXiv:astro-ph/0411220i. 

[17] G. L Rigopoulos and E. P. S. SheUard, Phys. Rev. D 68, 123518 (2003) arXiv:astro-ph/0306620j. 

[18] D. Langlois and F. Vernizzi, Phys. Rev. Lett. 95, 091303 (2005) |ar Xiv:astro-ph/0503416j . 

[19] D. Langlois and F. Vernizzi, Class. Quant. Grav. to appear (2010) 

[20] K. Koyama, arXiv:1002.0600" [hep-th]. 

[21] V. F. Mukhanov, Sov. Phys. JETP 67, 1297 (1988) [Zh. Eksp. Teor. Fiz. 94N7, 1 (1988)]. 
[22] M. Sasaki, Prog. Theor. Phys. 76, 1036 (1986). 



[23] C. Armendariz-Picon, T. Damour and V. F. Mukhanov, Phys. Lett. B 458, 209 (1999) (arXiv:hep- 
th/9904075|. 

[24] M. Alishahiha, E. Silverstein and D. Tong, Phys. Rev. D 70, 123505 (2004) |arXiv:hep-th/0404084|. 



[25] I. ZabaUa, Y. Rodriguez and D. H. Lyth, JGAP 0606, 013 (2006) |arXiv:astro-ph/0603534 



[26] C. T. Byrnes, K. Koyama, M. Sasaki and D. Wands, JCAP 0711, 027 (2007) [arXiv: 0705. 4096 

[hep-th]]. 

[27] R. Scoccimarro and J. Frieman, Astrophys. J. SuppL 105, 37 (1996) [arXiv:astro-ph/9509d47| . 



[28] M. Crocce and R. Scoccimarro, Phys. Rev. D 73, 063519 (2006) |arXiv:astro-ph/0509418 
[29] D. Seery, Class. Quant. Grav. to appear (2010). 



[30] L. Boubekeur and D. H. Lyth, Phys. Rev. D 73, 021301 (2006) |arXiv:astro-ph/0504046 



[31] C. T. Byrnes, M. Sasaki and D. Wands, Phys. Rev. D 74 (2006) 123519 [arXiv:astro-p h/0611075|. 

[32] C. T. Byrnes, S. Nurmi, G. Tasinato and D. Wands, arXiv:0911.278l) [astro-ph.CO]. 

[33] C. T. Byrnes, K. Y. Choi and L. M. H. Hall, JCAP 0902, 017 (2009) arXiv:0812.0807 [astro-ph]]. 

[34] M. Kawasaki, K. Nakayama, T. Sekiguchi, T. Suyama and F. Takahashi, JCAP 0811, 019 (2008) 



arXiv:0808.0009 [astro-ph]]. 



[35] D. Langlois, F. Vernizzi and D. Wands, JCAP 0812, 004 (2008) [arXiv:0809.4646 [astro-ph]] . 
[36] C. Hikage, K. Koyama, T. Matsubara, T. Takahashi and M. Yamaguchi, Mon. Not. Roy. Astron. 

Soc. 398, 2188 (2009) |arXi v:0812.3500l [astro-ph]]. 
[37] D. Langlois, Phys. Rev. D 59, 123512 (1999) |;aFXiv:astro-ph/9906080'. 
[38] F. Arroja and M. Sasaki, arXiv: 1002. 1376 [astro-ph.CO]. 



[39] D. H. Lyth and D. Wands, Phys. Rev. D 68, 103515 (2003) |arXiv:astro-ph/0306498 



[40] T. Falk, R. Rangarajan and M. Srednicki, Astrophys. J. 403, LI (1993) arX iv: astro-ph/920800 1 
[41] A. Gangui, F. Lucchin, S. Matarrese and S. MoUerach, Astrophys. J. 430, 447 (1994) |arXiv:astro- 
ph/93120331. 

[42] L. E. Allen, S. Gupta and D. Wands, JCAP 0601, 006 (2006) |arXiv:astro-ph/0509719|. 



[43] P. Creminelh and M. Zaldarriaga, JCAP 0410, 006 (2004) [arXiv:astro-ph/0407059 
[44] S. Mollerach, Phys. Rev. D 42, 313 (1990). 

[45] A. D. Linde and V. F. Mukhanov, Phys. Rev. D 56, 535 (1997) |arXiv:astro-ph/9610219|. 



[46] K. Enqvist and M. S. Sloth, Nucl. Phys. B 626, 395 (2002) |arXiv:hep-ph/0109214 



[47] D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) | arXiv:hep-ph/01 10002 . 
[48] T. Moroi and T. Takahashi, Phys. Lett. B 522, 215 (2001) [Erratum-ibid. B 539, 303 (2002)] 
arXiv:hep-ph/01 100961 . 

[49] D. H. Lyth, C. Ungarelli and D. Wands, Phys. Rev. D 67, 023503 (2003) |arXiv:astro-ph/0208055| . 



Local non- Gaussianity from inflation 



20 



B. A. Bassett, S. Tsujikawa and D. Wands, Rev. Mod. Phys. 78, 537 (2006) jarXiv:astro- 
ph/0507632|. 

M. Sasaki, J. Valiviita and D. Wands, Phys. Rev. D 74, 103003 (2006) |arXiv:astro-ph/0607627|. 



N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 69, 043503 (2004) |arXiv:hep-ph/0309033, . 

M. Giovannini, Phys. Rev. D 67, 123512 (2003) |arXiv:hep-ph/0301264\ 

K. Y. Choi and J. O. Gong, JCAP 0706, 007 (2007) arXiv:0704.2939 [astro-ph]]. 

H. AssaduUahi, J. VaUviita and D. Wands, Phys. Rev. D 76, 103003 (2007) [arXiv:0708.0223 

[hep-ph]]. 

Q. G. Huang, JCAP 0809, 017 (2008) |arXiv:0807.1567 [hep-th]]. 

K. Dimopoulos, G. Lazarides, D. Lyth and R. Ruiz de Austri, Phys. Rev. D 68, 123515 (2003) 

|arXiv:hep-ph/0308015] . 
K. Enqvist and S. Nurmi, JCAP 0510, 013 (2005) | arXiv:ast ro-ph/0508573] . 
Q. G. Huang, JCAP 0811, 005 (2008) |arXiv:0808.1793 [hep-th]]. 
K. Enqvist, S. Nurmi, O. Taanila and T. Takahashi, arXiv:0912.4657 [astro-ph.CO]. 



D. H. Lyth and D. Wands, Phys. Rev. D 68, 103516 (2003) arXiv:astro-ph/030650()] . 



S. Weinberg, Phys. Rev. D 67, 123504 (2003) ar Xiv:astro-ph/0302326 . 

E. Komatsu et ai, arXiv: 10 01. 4538 [astro-ph.CO]. 

D. Langlois and F. Vernizzi, Phys. Rev. D 70, 063522 (2004) [arXiv :astro-p h/0403258 V 

F. Ferrer, S. Rasanen and J. Vahviita, JCAP 0410, 010 (2004) [arXiv:astro-ph/0407300' . 

G. Lazarides, R. R. de Austri and R. Trotta, Phys. Rev. D 70, 123527 (2004) ^arXiv:hep- 
ph/0409335 |. 

K. Ichikawa, "T. Suyama, T. Takahashi and M. Yamaguchi, Phys. Rev. D 78, 023513 (2008) 
arXiv:080^T38 [astro-ph]]. 



M. Lemoine, J. Martin and J. Yokoyama, Phys. Rev. D 80, 123514 (2009) |a rXiv:0904.0126| [astro- 
ph.CO]]. 

G. I. Rigopoulos, E. P. S. Shellard and B. J. W. van Tent, Phys. Rev. D 73, 083522 (2006) 

arXiv:astro-ph/0506704 . 
L. Alabidi and D. H. Lyth, JCAP 0605, 016 (2006) |arXiv:astro-ph/051044l] . 



[41]. 
)9T 



F. Vernizzi and D. Wands, JCAP 0605 (2006) 019 [arXiv:astro-ph /0603799| 
D. Polarski and A. A. Starobinsky, Phys. Rev. D 50, 6123 (1994) [arXiv:astro-ph/9404061' . 
C. T. Byrnes, K. Y. Choi and L. M. H. HaU, JCAP 0810, 008 (2008) a rXiv:0807.110 1 [astro-ph]] 
T. Tanaka, T. Suyama and S. Yokoyama, arXiv:1003.5057 [astro-ph.CO]. 

C. T. Byrnes and K. Y. Choi, |arXiv:1002.3110( ]astro-ph.CO]. 

X. Chen, R. Easther and E. A. Lim, JCAP 0706, 023 (2007) [arXiv:astro-ph/0611645'. 

D. Langlois and L. Sorbo, JCAP 0908, 014 (2009) arXiv:0906.1813 [astro-ph.CO]]. 

F. Bernardeau, L. Kofman and J. P. Uzan, Phys. Rev. D 70, 083004 (2004) jarXiv:astro- 
ph/0403315|. 

D. H. Lyth, JCAP 0511, 006 (2005) 'arXiv:astro-ph/05 104431 . 
M. P. Salem, Phys. Rev. D 72, 123516 (2005) arXiv:astro-ph/0511146 . 
M. Sasaki, Prog. Theor. Phys. 120, 159 (2008) arXiv:0805.0974 [astro-ph]]. 
A. D. Linde, Phys. Rev. D 49, 748 (1994) [arXiv:astro-ph/9307002|. 
D. H. Lyth and A. Riotto, Phys. Rev. Lett. 97, 121301 (2006) arXiv:a stro-ph/0 607326j. 



G. Dvah, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 69, 023505 (2004) |arXiv:astro- 

ph/0303591,|. 
L. Kofman, arXiv:astro-ph/0303614' 

M. Zaldarriaga, Phys. Rev. D 69, 043508 (20 04) |arXiv:astro-ph/Q306006| . 



F. Vernizzi, Phys. Rev. D 69, 083526 (2004) [arXiv:astro-ph/0311167 



K. Kohri, D. H. Lyth and C. A. Valenzuela-Toledo, JCAP 1002, 023 (2010) [arXiv:0904.0793 

[hep-ph]]. 



G. Dvah, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 69, 083505 (2004) jarXiv:astro- 
|ph/0305548j . 



Local non- Gaussianity from inflation 



21 



[90] E. W. Kolb, A. Riotto and A. Vallinotto, Phys. Rev. D 71, 043513 (2005) |arXiv:astro- ph/0410546|. 
[91] L. Ackerman, C. W. Bauer, M. L. Graesser and M. B. Wise, Phys. Lett. B 611, 53 (2005) 
arXiv:astro-ph/0412007 . 

[92] E. W. Kolb, A. Riotto and A. Vallinotto, Phys. Rev. D 73, 023522 (2006) |arXiv:astro-ph/0511l"98]. 



[93] C. T. Byrnes and D. Wands, Phys. Rev. D 73, 063509 (2006) arXi v:astro-ph/0512195 

[94] T. Matsuda, JCAP 0703, 003 (2007) arXiv:hep-th/0610232|. 

[95] C. T. Byrnes, JCAP 0901, Oil (2009) |arXiv:0810.3913 [astro-ph] 

[96] M. Gasperini and G. Veneziano, Mod. Phys. Lett. A 8, 3701 (1993) [arXiv:hep-th /9309023 

[97] D. Wands, Phys. Rev. D 60, 023507 (1999) |arXiv:gr-qc/9809062 



[98] F. Finelh and R. Brandenberger, Phys. Rev. D 65, 103522 (2002) arXiv:hep-th/01 12249 



[99] Y. F. Cai, W. Xue, R. Brandenberger and X. Zhang, JCAP 0905, Oil (2009) [arXiv:0903.063T 
[astro-ph. CO]]. 

[100] J. Khoury, B. A. Ovrut, P. J. Stcinhardt and N. Turok, Phys. Rev. D 64, 123522 (2001) 
arXiv:h ep-th/0 103239| . 

R. Kallosh, L. Kofman and A. D. Linde, Phys. Rev. D 64, 123523 (2001) |arXiv:hep-th /0104073| . 
J. L. Lehners, P. McFadden, N. Turok and P. J. Steinhardt, Phys. Rev. D 76^ 103501 (2007) 

iarXiv:h ep-th/07 02153] . 
P. Creminelh and L. Senatore, JCAP 0711, 010 (2007) (arXiv:h ep-th/0702165| . 
D. H. Lyth, Phys. Lett. B 524, 1 (2002) arXiv:hep-ph/0106153|. 
A. Notari and A. Riotto, Nucl. Phys. B 644, 371 (2002) |arXi v:hep-th/0205019| . 
F. Finelh, Phys. Lett. B 545, 1 (2002) |arXiv:hep-th/0206112i 



[101 
[102 

[103 
[104 
[105 
[106 
[107' 

[108 
[109 

[110 

[111 
[112 

[113 
[114' 

[115 

[116 
[117" 

[118 
[119 

[120 
[121 
[122 
[123 



E. 1. Buchbinder, J. Khoury and B. ArOvrut," Phys. Rev. D 76, 123503 (2007) |arXiv:hep- 

th70702r54|. 

K. Koyama and D. Wands, JCAP 0704, 008 (2007) [arXiv:hep-th/0703040l 



K. Koyama, S. Mizuno and D. Wands, Class. Quant. Grav. 24, 3919 (2007) |arXiv:0704.1152 

[hcp-th]]. 



K. Koyama, S. Mizuno, F. Vernizzi and D. Wands, JCAP 0711, 024 (2007) [arXiv:0708.4321 

[hep-th]]. 

J. L. Lehners and S. Renaux-Petel, Phys. Rev. D 80, 063503 (2009) | arXiv:0906.053 0| [hep-th]]. 
E. I. Buchbinder, J. Khoury and B. A. Ovrut, Phys. Rev. Lett. 100, 171302 (2008) 

arXiv:0710.5172 [hep-th]]. 
T. Battefeld, Phys. Rev. D 77, 063503 (2008) |arXiv:0710.2540| [hep-th]]. 

A. Slosar, C. Hirata, U. Seljak, S. Ho and N. Padmanabhan, JCAP 0808, 031 (2008) 

arXiv:0805.3580 [astro-ph]]. 
J. Smidt, A. Amblard, A. Cooray, A. Heavens, D. Munshi and P. Serra, arXiv:1001.5026 [astro- 
ph.CO]. 

D. Wands and A. Slosar, Phys. Rev. D 79, 123507 (2009) arXiv:0902.W8l [astro-ph.CO]]. 

N. Dalai, O. Dore, D. Huterer and A. Shirokov, Phys. Rev. D 77, 123514 (2008) ^arXiv:0710.4560| 

[astro-ph]]. 

E. Komatsu et aZ., |arXiv:0902.4759| [astro-ph.CO]. 

N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Phys. Rept. 402, 103 (2004) |arXiv:astr(> 
ph/0406398 ]. 

C. Pitrou, J. P. Uzan and F. Bernardeau, |arXiv:1003.0481| [astro-ph.CO]. 
A. R. Liddle and D. H. Lyth, 

H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1 (1984). 
J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 



